Microscopic Derivation of Magnetic Coupling in Ca3Co206 
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Abstract 

For cobalt atoms placed along chain structures in Ca3Co2 06, we investigate the spin exchange coupling between 
atoms in high spin states. Consistent with experimental findings, the coupling is weakly ferromagnetic. 
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Quasi one-dimensional structures have long been 
acknowledged to exhibit intriguing transport and 
magnetic properties. New members in this fam- 
ily are Ca3Co20g [1] and a series of related iso- 
structural compounds, which crystallize in the i?3c 
structure [2]. In this structure, two inequivalent 
cobalt atoms exist, one in an octahedral environ- 
ment (labeled Col) , and the other one in a trigonal- 
prismatic environment (labeled Co2). Both atoms 
are in a 3+ configuration, but in different spin 
states. Cobalt atoms form covalent bonds along 
chains while the interchain coupling is far weaker. 
One therefore expects that this rather unique sit- 
uation is going to give rise to a peculiar magnetic 
response. 

Along the chains the cobalt spins tend to order 
ferromagnetically whereas the interchain coupling 
is antiferromagnetic. However, the antiferromag- 
netic order is affected by the frustration inherent 
in the crystal structure. At low temperature, this 
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system appears to be in a partially disordered an- 
tiferromagnetic state. On top, the magnetization 
versus field curves display several plateaus [1] . 

An appealing starting point to the understand- 
ing of the magnetic structure of a particular com- 
pound is provided by the Kanamori-Goodenough- 
Anderson rules [3]. Given the electronic configu- 
ration of the ions one can estimate their magnetic 
couplings and, at mean-field level, the magnetic 
ground state. In Ca3Co2 0g, the problem is more 
involved, in particular since each second Co 3+ is 
non-magnetic, as indicated by neutron scattering 
data [2]. One therefore needs to determine the 
next-nearest neighbor magnetic coupling. This is 
the purpose of the paper. 

To that aim one needs i) to determine the elec- 
tronic configuration of the ions, and ii) the ma- 
trix elements of the hopping operator. Both are de- 
scribed by a set of parameters obtained from den- 
sity functional theory calculations using the aug- 
mented spherical wave method ( ASW) [5] . In par- 
ticular, it turns out that both Co ions are in a 3d 6 
configuration. The parameters entering the local 
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part of the standard multi-band Hubbard Hamil- 
tonian, upon integrating out the oxygen degrees 
of freedom, are the Hubbard U, the Hund's rule 
coupling Jh, and the crystal field splitting 10 Dq. 
For Co2, the splitting of the "three-fold degener- 
ate" lowest level turns out to be small [6], and is 
neglected. In contrast to 10 Dq both U and Jh are 
weakly affected by the difference between the pris- 
matic and octahedral environments. The six elec- 
trons can be in three spin configurations: a low 
sp in one (LS) (S = 0) wit h energy E LS = (30f7 - 
y/(3J H + 20Dq) 2 + 2W H - 57 J H - 28Dq)/2, an 
intermediate spin one (IS) (5 = 1) with energy 
Eis = (SOU - ^{J H + 20Dq) 2 + 8J H - 65 J H - 
8Dq)/2, or a high spin one (HS) (S = 2) with en- 
ergy E HS = 15U— 38 J H - ^Dq. As a result the LS 
configuration is stable for Jh < 2.68Dq, and the 
HS one otherwise. Stabilizing the IS configuration 
requires different physical input, such as an addi- 
tional splitting of the e g level or a lowering of the 
symmetry [4]. 

We are now in the position to estimate the 
intra-chain magnetic couplings. For Col we choose 
the coordinate axes to point towards the oxygens, 
while for Co2 we choose the z-axis to point along 
the c-direction of the i?3c group. Therefore when 
determining the matrix elements of the hopping 
term one needs to carefully distinguish the rele- 
vant contributions. It turns out that the 3z 2 — r 2 
orbital on Co2 identically couples to the x'y', x'z', 
and y'z' orbitals on Col, with amplitude t, while 
other matrix elements are negligible. At this stage 
one can determine the magnetic ground state per- 
forming a hopping expansion. It is obvious that 
the magnetic structure does not enter the result to 
order t 2 . To order t 4 the exchange paths between 
two Co2 sites are sensitive to the magnetic struc- 
ture and the respective Hund's rule coupling on 
Col. A typical exchange path is sketched in Fig. 1 
for a ferromagnetic configuration. Any tig down 
spin on Col can hop onto the 3z 2 — r 2 orbital on 
one neighboring Co2, and any remaining ti g down 
spin on Col can hop onto the 3z 2 — r 2 orbital onto 
the second neighboring Co2, (this leaves Col on a 
spin state where Jh is important), and then both 
electrons return to Col. There are 48 such contri- 
butions. For an antiferromagnetic configuration, 
we find that: i) there are only 42 such paths, and 
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Fig. 1. Exchange paths, as discussed in the text. 

ii) it leaves Col on an intermediate low spin state, 
which has a higher energy. As a result the ferro- 
magnetic configuration has a lower energy. Using 
the numerical values resulting from the ASW 
calculation [6] and [7 = 5.5 eV, we finally obtain 

Ep - E AF ~ ~10- 3 t 4 /eV 3 . (1) 

Using t ~ 0.8 eV, the energy gain is around 40 K, 
in good agreement with the onset temperature of 
the ferromagnetic ordering along the chains. 

In summary, we have identified a mechanism 
leading to a ferromagnetic coupling between mag- 
netic ions through a non-magnetic one. The driv- 
ing forces are i) the Hund's rule coupling on the 
non-magnetic ion, and ii) the large number of con- 
tributions to the magnetic energy due to the R3c 
structure. The magnetic coupling is found to be in 
good agreement with experiment. 
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